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Near-Field Noise Computation for a Subsonic Coannular Jet 


Ching Y. Loh 
Taitech, Inc. 

Beavercreek, Ohio 45430 

Lennart S. Hultgren and Philip C.E. Jorgenson 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

A high-Reynolds-number, subsonic coannular jet is simulated, using a three-dimensional finite-volume LES 
method, with emphasis on the near field noise. The nozzle geometry used is the NASA Glenn 3BB baseline 
model. The numerical results are generally in good agreement with existing experimental findings. 

I. Introduction 

The reduction of aircraft noise is an important issue due to environmental regulations and economic consequences 
of airport community noise. For modern commercial aircraft engines at subsonic speeds, turbulent mixing noise 
generated in the jet flow is, in general, a major noise component and is dominant at takeoff. The study of jet noise 
has been a vital part of aeroacoustics for over half a century and continues to be a topic of many experimental and 
theoretical investigations; Refs. 1-3 provide a comprehensive discussion and further references. A rapid advance in 
its understanding has occurred during the past decade, perhaps in conjunction with the development of computational 
aeroacoustics (CAA) and of improved measurement techniques. The issues, challenges, and contributions of CAA are 
discussed in Refs. 4-6 and Refs. 7-12 provide examples of advances in experimental methods. 

Fine scale turbulence in jet shear flow produces sound propagating predominantly in the lateral (90°) direction, 13-1 6 
whereas the convection/evolution of large-scale turbulence structures in the jet produces sound propagating down- 
stream at angles close to the jet axis. The work in this paper concerns the sound produced by the large-scale structures 
in subsonic jets. 

Substantial progress in numerical computation of jet mixing noise associated with large-scale structures has been 
made during recent years by using large-eddy simulation (LES) techniques. Mankbadi et al. 17 performed an axisyrn- 
metric LES computation of a supersonic, high-Reynolds-number jet and determined its far-field radiated sound using 
an acoustic analogy technique. Zhao et al., 18 Bogey et al., 19 and Uzun et al. 20 computed jet flow and the associated 
noise by using fully three-dimensional LES at a Mach number of 0.9. The Reynolds numbers in these simulations were 
low to moderate; 5, 000, 65, 000, and 100, 000, respectively. DeBonis and Scott 21 and Loh and Hultgren 22 used LES 
to calculate the hydrodynamics of a high-Reynolds-number, supersonic jet of Mach number 1.4. 8,9 More recent work 
and discussions of the state-of-the-art can be found in Refs. 23-26. Except for our own work, 22 which utilized the 
second-order (in space and time) finite-volume CE/SE Navier-Stokes solver, the numerical schemes used by the above 
were of high-order finite-difference type. Recent LES studies of coaxial jets have also been carried out by Refs. 27,28. 

The present paper utilizes an enhanced time-accurate, upwind, finite-volume scheme for compressible flows on 
unstructured grids. 29 The method (ETAU) is nominally second-order accurate in space and time and is capable of 
yielding high resolution in the presence of sharp gradients. A built-in dissipation model helps to overcome potential 
pathological behaviors 30,31 of upwind schemes and renders the ETAU scheme robust and viable for practical compu- 
tations. 


II. The Jet Flow 

The current computation is for a subsonic dual-stream flow issuing from the NASA Glenn 3BB baseline coan- 
nular nozzle, 32 see Fig. 1, into a co-flow stream of Mach number 0.28. The nozzle operational conditions corre- 
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spond to the ones used by Koch et a l? 2 in their comparison of numerical RANS results with experimental data. 

In what follows, the exit diameter D of the of the core 
nozzle, the density p^, the speed of sound a^, and 
the (static) temperature T x in the co-flow are used as 
scales to make all the dependent variables nondimen- 
sional. Ideal nozzle flow is fully determined by spec- 
ifying the back-pressure ratio (ambient static pressure 
divided by stream total pressure) and the total temper- 
ature of the stream. To determine the plenum conditions, 
the only additional information needed is the ratio of the 
plenum and exit areas. 

Using the ideal gas isentropic relations it follows that 
for non-choked flow in a nozzle 

1 < Ptot/Poo < Ptot/P * = [(7 + l)/2] 7/(T_1) (1) 


Figure 1. NASA Glenn baseline nozzle, 3BB 12 

Pe = Poo = 1/7) Pe = 1 /T e , T e 


where P to t, Poo, and P* denote the total, back, and sonic 
pressures, and 7 = C p /C v is the ratio of specific heats. 
Within the restrictions of Eq. 1 , it follows that the nondi- 
mensional nozzle exit conditions are 

= Ttot(P 00 /Pt 0 t) {r> - 1)h , |u e | = M e s/r e , (2) 



where p, T, and u = (u, v, w) denote density, temperature, and velocity vector; and the exit Mach number is given by 


Me = y^pWPoo^-U/T _ i ]/( 7 - 1 ) . 

Conservation of mass in combination with the ideal gas isentropic relations can be used to show that 


7>+l 

Mr/ [1 + 1(7 - 1 )M, 2 ] 


~l~ 1 

(A e /A r )M e / [1 + 1( 7 - 1 )M e 2 ] ^ ; 


( 3 ) 

( 4 ) 


where A r /A e is the ratio of the plenum (reservoir) and exit areas. Since M e is known from Eq. 3 for a given back- 
pressure ratio, the right-hand side of Eq. 4 is also known, once the area ratio is specified, and hence the plenum Mach 
number, M r can be determined iteratively. The nondimensional plenum conditions are then obtained using 

Pr = Ptot/[l + \(l-l)M /:]^ 1 , p r = 1 P r /T r , T r = T tot / [l + | (7 — 1)M/] , Kl =M r y/r r . (5) 


Eqs. 3-5 together with the nozzle conditions in Table 1 fully define the mean flow. 


Table 1. Nondimensional Nozzle Flow Conditions 


Core Stream 

Fan Stream 

Free Stream 

Ptot/Poo 

1.6801 

Ptot/Poo 

1.8301 

Ptot./ Poo 

1.0560 

Ttot/To o 

2.8320 

Ttot /Too 

1.1328 

Ttot /Too 

1.0157 

A r /A e 

2.7762 

A r / A e 

2.9327 

Moo 

0.28 


III. Numerical Method 

The numerical scheme adopted in this work is an extension to unstructured hexahedral grids of the low-dissipation, 
three-dimensional, finite-volume ETAU scheme. 29 It is applied to an LES version of the compressible Navier-Stokes 
equations. The ETAU scheme has formal second-order accuracy in both space and time. It has an additional small 
amount of omni-directional numerical dissipation, compared to standard (see Ref. 33) upwind schemes, that renders it 
capable of handling a broad range of flow regimes. 
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A. Computational Domain and Unstructured Hexahedral Grid 



Figure 2. Geometry of the computational domain; the core-nozzle exit diameter is used as length scale. 


The computational domain including the coannular 
nozzle is sketched in Fig. 2, with a close-up in Fig 3. It 
ranges from x = —3.64 to x = 18.20, where x is the 
streamwise coordinate. The domain diameter is about 
11. The numerical simulation begins at the inlets of the 
coannular nozzle. The fan nozzle exit is at x = —1.06; 
the core nozzle exit is at x = 0; and the centerbody ends 
at x = 0.97. 

For fully three-dimensional computations of compli- 
cated problems such as the current jet noise problem 
where several distinct and disparate length scales need to 
be resolved, 4 a large number of computational cells are 
needed to provide sufficient numerical resolution. The 
choice of an unstructured hexahedral grid for the ETAU 
method leads to a significant reduction of the number 
of computational cells, with an associated reduction in 
computer memory and CPU time requirements, as com- 
pared to the commonly used tetrahedral one. The draw- 
back is a modest increase in numerical dissipation. 

The unstructured hexahedral grid is generated using 
the Gridgen 34 software package. The streamwise step 
size Ax ranges from 0.022 at the nozzle exit (x = 0) to 
about 0.05 at x = 16.2. Between x = 16.2 and x = 18.2 there is a buffer zone in which Ax increases exponentially 
in size to provide high numerical damping. In total there are about 5.36 million cells in the computational domain, 
including the buffer zone. The domain covers the near-field region which is the source region for the noise. 

Employing a hexahedral grid instead of a tetrahedral one for the current ETAU scheme helps to reduce the number 
of cells from tens of millions to a few millions. Still, the number of computational cells is very large and the simulation 
is implemented using parallel-computation techniques. To accomplish this, the unstructured hexahedral grid generated 
for the entire computational domain is decomposed into subdomains (128 at present) using the efficient mesh parti- 



Figure 3. Close-up of Computational Grid Plane 
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tioning code METIS. 35 The ETAU flow solver uses the message-passing-interface (MPI) library to exchange pertinent 
information between neighboring subdomains during the simulation. 

B. Conservation Form of the Three-Dimensional Unsteady Compressible Navier-Stokes Equations 

The Navier-Stokes equations for a perfect gas can be written in the following nondimensional vector form, e.g. Ref. 36: 

U t + F x + G y + H z = 0, (6) 

where x, y, and z are the streamwise, transverse, and spanwise coordinates, and t is time. The five components of the 
conservative flow variable vector U are the density, the three momentum components, and the total energy per unit 
volume, i.e., 

Ui = p, U 2 = pu , Uz = pv, U4 = pw , 

U 5 = p/(l ~ 1) + piu 2 + v 2 + w 2 )/2, 

where, as before, p, u, and p denote the density, streamwise velocity, and static pressure, and v and w denote the 
transverse and spanwise velocities. The flux vectors in the x, y, and z directions, F, G, and H , respectively, have 
both inviscid and viscous contributions and the details of these flux components may be found in a textbook 36 or in a 
previous paper 37 by the present authors. 

The nondimensional equations above could form the basis for a direct numerical simulation (DNS). Currently, 
this is not feasible for the high-Reynolds-number flows of interest here. The present computations are of the LES 
type, where the effects of the unresolved scales, i.e., scales smaller than the grid resolution, are modeled. In this case, 
the equations above are considered the filtered ones governing the resolved scales, where p and p are interpreted as 
simply averaged and the other dependent variables are interpreted as Favre, (i.e., density-weighted) averaged over the 
subgrid scale. The simplest model for the effects of the unresolved scales, or subgrid scales (SGS), on the motion is 
obtained by using a Boussinesq eddy-diffusivity assumption for the subgrid shear stresses and heat flux coupled by 
Smargorinsky’s model for the eddy viscosity (eg. Ref. 38). The nondimensional viscosity p, occurring in the viscous 
contribution to the flux terms in Eq. 6, is then replaced by 

p= 1/Re +(C s S) 2 (S ij S ij ) 1 / 2 , (7) 

where S l3 = \ is the Favre averaged rate-of-strain tensor, S = ( 5x6yd z ) 3 / 3 is a local measure of the 

grid size, C s is a (nondimensional) constant, and Re = CLooD/v <*, is the Reynolds number for the computation. The 
inverse Reynolds number 1 / Re = 1.8 x ICC 6 in the present computations. Also, note that the simplifying assumption 
is also made that the subgrid-scale Prandtl number equals the laminar one, which is taken to have the value 0.72. 

Note that Fureby and Grinstein 39 point out that applying flux limiters to finite-volume methods (as is done in the 
current ETAU scheme), even in the absence of any explicit LES assumptions, effectively leads to LES schemes with 
minimal implicit SGS models. They demonstrated through an ‘error’ analysis of a particular scheme that the flux- 
limiters (essentially low-pass frequency filters) build into the algorithm produces additional terms in the equivalent 
differential forms of the momentum and energy equations that can be interpreted as the SGS stress tensor and flux, 
respectively. Hence, there is an additional implicit SGS model inherent in the particular ETAU scheme used in addition 
to the explicit assumption in Eq. 7. 

C. Enhanced Time-Accurate Upwind Scheme 

By considering (x, y, z, t ) as coordinates of a four-dimensional Euclidean space, E±, and using the Gauss divergence 
theorem, it follows that Eq. 6 is equivalent to the following integral conservation law: 

(f I m dS = 0, m= 1,2,..., 5, (8) 

JS(V) 

where S(V) denotes the surface around a volume V in E± and I rn = (F m , G m , H m , U m ). These equations are 
discretized, using an unstructured hexahedral grid, with the flow variables stored at each cell center. 

Consider the flow variable U and its spatial derivatives, U x , U v , and U z , as known at the cell centers at time level 
n. The evolution step of the enhanced time-accurate upwind scheme then advances the value at a cell center i to time 
level n + 1 through 

U\ n+1) = /3U Z + (1 - f3)U- n) - f>!" +1/2)r U* + G^ +1/2) n Vk + H^ +1/2) n zk ]As k , (9) 

v fc= 1 
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where Fk, Gk, Hk, n xk , n Vk , n zk , and A s k , are the three flux vectors evaluated for half a time step at the centers of 
the twelve triangular surfaces that define the hexahedral cell, the corresponding out-going unit normal vector compo- 
nents, and areas, respectively; Av is the volume of the hexahedral cell and At is the time step; Ui is an average based 
on the values in the six surrounding cells at time step n and f3 is a parameter (0 < /3 < 1). Here, U, is taken to be the 
simple average of the 12 right states at time level n, 

Ui = ^f^(ul n) +U^Ax k + U^Ay k + U^Az k ) j = j(k), (10) 

fc= l 

where the right states are linear Taylor series extrapolations from the corresponding neighboring cell center, j(k), to 
the surface center, k. For /3 = 0, Eq. 9 reduces to the usual update formula for an upwind scheme that provides the 
cell-averaged values at the next time level. The flux vectors are generated using a Riemann solver, where the required 
right and left states are obtained by linear Taylor expansion, in both space and time, and Cauchy-Kowaleswki/Lax- 
Wendroff time stepping. 29 Each right state involves only the appropriate neighboring cell and the corresponding left 
state only involves the current cell. The excluded case /? = 1, in a sense, corresponds to the previously used 22 ' 40 
CE/SE scheme, except of course a Riemann solver is not used to obtain the surface flux vectors in that method. The 
value /3 = 0.01 is used throughout the present work. 

The reconstruction step of the ETAU scheme then obtains the flow variable U as a linear vector function within 
each cell i at the new time level n + 1 by determining the spatial derivatives at each cell center, U ( :) n ' +1 \ U^ +l \ and 

U i" +1) • This procedure uses the updated flow variable information at the current cell center and its six adjacent cells 
and leads to an over-determined system. Similarly to Ref 29, a multidimensional van Albada limiter 41 is applied to 
obtain the required solution. The evolution/reconstruction procedure can then be repeated to march the solution further 
in time. 


D. Initial and Boundary Conditions 

Initially, the flow of the entire computational domain is set at the ambient flow conditions, i.e., 

Po = 1 , Po — 

7 

u 0 = Moo, V 0 =0, w 0 = 0. 

At the exterior upstream inlet boundary, the conservative flow variables are specified to be the same as the ambient 
flow and their spatial derivatives are held at zero. At the nozzle inlets, the reservoir conditions as described in § II are 
imposed. Of course, v r ,w r , and all spatial derivatives are held at zero. 

Non-reflecting boundary conditions 42,43 are imposed at the circumferential and outflow boundaries. The slip 
boundary condition is applied on all the nozzle walls, external as well as internal, except on the external centerbodies 
for both fan and core streams where the no-slip condition is applied. 


E. Artificial Forcing 

As is common, e.g. Ref 19, in numerical studies of jet mixing noise, the simulation allows for a low-level, artificial 
excitation, or turbulence seeding, of the initial jet shear layers. This seeding is accomplished by adding divergence- 
free velocity perturbations inside the jet shear layers near the nozzle lips. The disturbances are temporally random and 
are of the vortex-ring form used in, 19 namely 


ll 


/ 



T(x,r,0\ x p , 


r P ), 


( 11 ) 


X-X p , / cos e 

' smf 


( 12 ) 


where 

25 r 16 

F{x,r,9;x p ,r p ) = " £ p exp (-rf 2 ln2/A 2 ) ^ e n cos(n9 + (j> n ) (13) 

A P n=l 

and r = ( y 2 + z 2 ) 3 / 2 and 9 = sin _1 (;z/r) are the radial and azimuthal coordinates; x p , r p and A p = 0.02 define 
the location and spatial extent of the perturbation, and d = [{x — x p ) 2 + (r — rp) 2 ] 1 / 2 ; e n and <j) n are uniformly 
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distributed random numbers between —1 and 1, and 0 and 2tt, respectively. They represent the random perturbation 
amplitude and phase, and are independently generated for the two shear-layers at every time step. For the co-flow/fan 
shear layer, x p = —0.99 and r p = 1.0; and for the fan/core shear layer, x p = 0.05 and r p = 0.5. The overall 
perturbation amplitude S p is, in general, set to 0.003 in the present computations. This value is representative of 
a typical disturbance environment for laminar-nozzle flow. The perturbations above are then directly added to the 
velocity field. 



Figure 4. Instantaneous distributions in a spanwise plane at t = 440 , random forcing; top, streamwise velocity; bottom, spanwise vorticity 


IV. Numerical Results 

The simulation starts from the uniform initial conditions described above with the artificial forcing turned on using 
an amplitude S p = 0.003. The simulation is first run to t = 100, or 100, 000 time steps ( At = 0.001), to ensure that all 
start-up transients are convected out of the computational domain and that the jet flow is fully established. This time, 
t = 100, corresponds roughly to five acoustic times. The acoustic time is a measure of how long it takes a sound wave 
to pass through the external computational domain. The simulation is then run for an additional 340, 000 time steps 
to t = 440, or approximately an additional 17 acoustic times, to achieve sufficient accuracy in the spectral analysis of 
time series data that is saved at various locations in the computation domain every 20 time steps. 
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A. Instantaneous Flow Fields 


Figure 4 shows the computed instantaneous streamwise velocity (top) and spanwise vorticity (bottom) fields in the 
spanwise plane, z = 0, at t = 440. The instantaneous vorticity field, in particular, indicates that the outer core starts to 
disappear, or merge with the inner core at about 5 diameters downstream of the nozzle exit. The inner potential core 
ends around 8 diameters downstream of the nozzle exit and the mixing process then rapidly gains strength. 



Figure 5. Instantaneous isobars at t = 440, random forcing 


Figure 5 shows the computed instantaneous three-dimensional pressure isosurfaces at t = 440 for two levels, 
p = 0.71 and p = 0.715; the nondimensional ambient pressure is 1/1.4. The isosurfaces show downstream radiating 
three-dimensional sound waves caused by the large-scale coherent structures in the jet flow, i.e., eddy-Mach wave 
radiation. Note that even though both the core and fan jets are subsonic, with theoretical Mach numbers of 0.90 and 
0.97, their theoretical acoustic Mach numbers are 1.4 and 0.95. The computations show that the unsteady motion is a 
truly three-dimensional flow. 

B. Time-Averaged Data 

Figure 6 shows the radial variation of the time-averaged streamwise velocity for five streamwise stations, x = 1.11, 
5.12, 8.42, 11.73, and 15.03. The solid lines and symbols represent the numerical results and the experimental data 
presented in Ref. 32, respectively. The computed result is in very good agreement with the experimental data at the first 
streamwise station indicating that the early jet shear layers are quite well resolved in the simulation. The agreement is 
also very good at the fourth station. At this station, x = 11.73, the initial dual-core jet has evolved into a single-stream 
jet. Ref. 32 found experimentally that the jet flow is a self-similar single stream at locations greater than about six fan 
diameters downstream of the fan nozzle exit plane, which corresponds to about x = 11 here. The numerical results 
are in general agreement with the experimental results for stations two, three, and five. The slight asymmetry clearly 
visible in the distribution for the last station indicates that a longer observation time is needed to fully capture the mean 
flow at this location. It may also indicate that not enough energy in the large-scale motions is being cascaded down to 
the small-scale motions and dissipated. 

Figure 7 compares computed and experimental values for the radial distribution of rms streamwise velocity for 
three, x = 1.11, 11.73, and 15.03, of the streamwise stations above. The apparent discrepancy at the first station is 
because, at this location, the unsteady jet flow is still dominated by small scales, which of course can not be accounted 
for in the current LES simulation. The numerical values are representative of the applied forcing at the nozzle lip. 
However, the agreement between simulation and experimental results is good at the last two streamwise stations 
where the unsteady jet flow is dominated by large-scale structures. The asymmetry in the numerical result, that is also 
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Figure 6. Radial variation of time-averaged streamwise velocity at five streamwise locations; symbols: experiment, 32 solid lines: LES 



Figure 7. Radial variation of rms streamwise velocity at three streamwise locations; symbols: experiment, 32 solid lines: LES 


evident in the computed mean-flow profile above, is further evidence that the observation time at the last station is not 
long enough. 


V. Concluding Remarks 

The present mixing-noise computation for a subsonic coannular jet demonstrates that the near-field large-scale 
unsteady flow can be determined using a low-order (2nd order) finite-volume LES scheme, such as the current one. 
Instantaneous iso-surfaces of pressure show the radiating sound waves originating from the large-scale structures in the 
jet. These as well as instantaneous streamwise-velocity results demonstrate the highly unsteady and three-dimensional 
nature of the flow. Generally, the time-averaged results, such as for the streamwise velocity and rms streamwise 
velocity profiles agree well with existing experimental data. 
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